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Abstract 

Using the steepest descent method for oscillatory Riemann-Hilbert problems introduced by Deift and 
Zhou [Ann. Math. 137(1993), 295-368], we derive asymptotic formulas for the Meixner polynomials in 
two regions of the complex plane separated by the boundary of a rectangle. The asymptotic formula 
on the boundary of the rectangle is obtained by taking limits from either inside or outside. Our results 
agree with the ones obtained earlier for z on the positive real line by using the steepest descent method 
for integrals [Constr. Approx. 14(1998), 113-150]. 
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1 Introduction 



In this paper, we investigate the asymptotic behavior of the Meixner polynomials. These polynomials have 
many applications in statistical physics. For instance, they are used in the study of the shape fluctuations 
in a certain two dimensional random growth model; see [10] and the references therein. 
For /3 > and < c < 1, the Meixner polynomials are given by 

m„(z;/3,c) = (/3)„-2Fi(-n,-^;/3;l-c-i). (1.1) 

They satisfy the discrete orthogonality relation 

f2mnik-l3,c)mpik;f3,c)^^ = (1 - c)-^c-"n!(/3)„ V (1.2) 

fe=0 

This notation is adopted in ^ § 10.24] and also in [8]. 

Using probabilistic arguments, Maejima and Van Assche |llj have given an asymptotic formula for 
m„(na; f3, c) when a < and is a positive integer. Their result is in terms of elementary functions. By 
using the steepest-descent method for integrals, Jin and Wong [8] have derived two infinite asymptotic 
expansions for m„(na; /?, c); one holds uniformly for 0<e<a;<l + e, and the other holds uniformly for 
1 — e<a<M<cxD. Both expansions involve the parabolic cylinder function and its derivative. 

In view of Gauss's contiguous relations for hypergeometric functions [1, § 15.2] and the connection 
formula [8] 

m„(-x - /?;/?, c^^) = c"-mn{x; c), 

we may restrict our study to the case 1 < /3 < 2 and < c < 1. Fixing < c < 1 and 1 < /3 < 2, 
we intend to investigate the large-n behavior of mn{nz — /3/2;/3,c) for z in the whole complex plane, 
including neighborhood of the origin and regions extending to infinity. Our approach is based on the 
steepest-descent method for oscillatory Riemann-Hilbert problems, first introduced by Deift and Zhou [6] 
for nonlinear partial differential equations, later developed in |5j for orthogonal polynomials with respect 
to exponential weights, and further extended in [2l|3] to a general class of discrete orthogonal polynomials. 

A direct application of the method in [21 13] would, however, only give local asymptotics. For instance, 
in the case of Meixner polynomials, one would have to divide the complex plane into at least six regions 
(one near the origin, two near the two turning points and three in between, including an unbounded one), 
and give correspondingly six different asymptotic formulas. To reduce the number of these regions, we 
shall make some modifications to the method in [2, 3]. Our approach is motivated by the previous work 
in [4, 13, 15, 17, 18], and the main idea is to extend, as large as possible, the two regions of validity 
of the two asymptotic formulas near the two turning points. There have already been several examples 
in which we only need two regions with appropriate asymptotic formulas to cover the entire plane; the 
Hermite polynomial [17] is one of such examples. However, for discrete orthogonal polynomials, there 
might be cuts starting from the finite endpoints of the intervals of orthogonality. For instance, in the case 
of Krawtchouk polynomials considered in [4], there are two cuts (— oo,0] and [l,oo) where no asymptotic 
formulas are given. In the present paper, we shall give two asymptotic formulas for the Meixner polynomial 
mn{nz — /3/2;/3,c), one valid inside a rectangle with two vertical boundary lines passimg through z = 
and z = 1, and the other valid outside the rectangle. Both formulas can be extended slightly beyond the 
boundary of the rectangle, and they are asymptotically equal to each other in the overlapping region. 
The material in this paper is arranged as follows. In Section 2, we use a standard method to relate the 
Meixner polynomials to a Riemann-Hilbert problem for a matrix-valued function. The motivation and 
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details of this standard procedure can be found in [21 [3] and the reference given there. In Section 3, we 
introduce some auxihary functions which will be used in Section 4 for the construction of our parametrix. 
In Section 4, we also prove that this parametrix is asymptotically equal to the solution of the Riemann- 
Hilbert problem formulated in Section 2. In Section 5, we state our main result and make the remark 
that our formulas agree with the ones already existing in the literature. 



2 Standard formulation of Riemann-Hilbert problem 

From (jl.ip . we note that the leading coefficient of m„(z;/3,c) is (1 — c^^)". Thus, the monic Meixner 
polynomials are given by 



7r„(z) := (1 - -) "m„(z;/3,c). 



(2.1) 



For convenience, in (2.1) we have suppressed the dependence of '7r„(z) on c and /3. Furthermore, throughout 
the paper we shall fix the parameters c S (0, 1) and /3 E [1, 2). The orthogonality property of 7r„(z) can 
be easily derived from (II. 2p . and we have 



^TTn(,k)7rp{k)w{k) = 5nphl, 



k=0 



where 



and 



7^ 



(1 - c)2"+/5c-" 

r(n + /3)r(n + l) 



r(z + /3) , 

wiz) := — -c . 

^ ' r(z + i) 



Let P{z) be the 2x2 matrix defined by 



Piz) 



'Piiiz) Puiz) 

Mz) P22{z) 



T^niz) 



g 7rn{k)w{k) \ 
k=o z — k 



,jl_lTTn-l{z) 

\ k=0 



7n-l^n-l(fc)w(fe) 
Z-k J 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



A proof of the following result can be found in [3^ Section 1.5.1]. The only difference is that their N 
should be replaced by oo. 



Proposition 2.1. The matrix-valued function P{z) defined in /12. 5\) is the unique solution of the following 
interpolation problem: 



(PI) Piz) is analytic in C \ N; 
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(P2) at each z = A; G N, the first column of P{z) is analytic and the second column of P{z) has a simple 
pole with residue 

(P3) for z bounded away from N, P{z) ( ^ „ j = / + Odzl^-*^) as z ^ oo. 
Let X denote the set defined by 

X:={X,}glo, where X, := (2.7) 
cf. [3], [4] and [13]. The X^s are called nodes. For the sake of simplicity, we put 

n-l 

B{z):=J{{z-X,). (2.8) 

3=0 

Our first transformation is given by 

Q{z) := n-"'^='P(nz-/3/2)5(z)-'^3 

where cr^ := is a Pauli matrix. In this paper, we shall also make use of another Pauli matrix, 

namely ai := ; see Section 4. 

Proposition 2.2. The matrix-valued function Q{z) defined in \2. 9\) is the unique solution of the following 
interpolation problem: 

(Ql) Q{z) is analytic in C \ X; 

(Q2) at each node with A; G N and k >n, the first column of Q{z) is analytic and the second column 
of Q{z) has a simple pole with residue 

z=Xk z^Xk yj U J 

at each node X^ with A; G N and k <n, the second column of Q{z) is analytic and the first column 
of Q{z) has a simple pole with residue 



( o\ 



Res Q[z^ = lim Q{z) 



2 



{z - Xk) 
\w{nz- (3/2)B{zf V 



(2.11) 
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(Q3) for z hounded away from X, Q{z) = I + 0{\z\ as z ^ oo. 

Proof. This is obvious from Proposition 12.11 and the definition of Q{z) in (j2.9p . ■ 

The purpose of our next transformation is to remove the poles in the interpolation problem for Q{z) 
(cf. [31 Section 4.2]). Let 5 > be a sufficiently small number. We define (see Figure [T] below) 

R{z):=Q{z)( ^\ (2.12a) 



-A±(z) 1 



for Rez G (0,1) and itlmz G (0,(5), and 



R{z) := Q{z) {I ^f'A (2.12b) 



for Rez G (l,oo) and iblmz G (0,5), and 



for Rez ^ [0, oo) or Imz ^ [—(5,15], where 



R{z) := Q{z) (2.12c) 



miw{nz- p/2)B{zf 
•= eTi-(n.-/3/2) sin(n7rz-/37r/2)' ^^'^^^ 
g±i7r(n^-/3/2) sin(n7rz - /37r/2) 
^^^^^ n™(nz-/3/2)i?(z)2 ' ^^"^^^ 

Lemma 2.3. For eac/i /c G N, i/ie singularity of R{z) at the node = ^"^^^^ is removable, that is, 
Res R{z) = 0. 

Proof. For any A; G N with k > n, we have = ^^^"^ > 1 since 1 < /3 < 2. From (2.13), it is evident 
that the residue of V± (z) at z = is 

Res V±(z) = w{nXk - ^/2)B{Xk f. 

From (|2.10p . we also note that the residue of Qi2{z) at z = Xj. is Qii(Xfc) multiplied by w{nXk — 
P/2)B{Xkf. Thus, it follows from (2.12b) that the residue of Ri2{z) = Qi2{z) -V ±{z)Qii{z) ai z = Xk 
is zero. Similarly, one can show that Res R22{z) = 0. Since Rii{z) = Qii{z) and i?2i(-z) = Q2i{z), and 

z=Xf, 

since Qii{z) and Q2i{z) are analytic by Proposition 2.2, the residues of Rii{z) and R2i{z) at are zero. 
For any A; G N with k < n, we have Xk = < 1 since 1 < /3 < 2. From ([23]), ([2TT]1 and (pTi]) . we 

observe that 

n— 1 

Res A±(z) = — — —-Y\{Xk-Xj)-^, 
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Thus, the residue of Rii{z) = Qu{z) — A±(z)Qi2{z) at z = Xj. is zero. Similarly, one can prove that the 
residue of R2i{z) at z = Xj^ is also zero. Since Ri2{z) = Qi2{z) and R22{z)=Q22{z) , and since Qi2{z) and 
Q22{z) are analytic by Proposition 2.2, the residues of Ri2{z) and R22{z) at X^ are zero. This completes 
the proof of the lemma. ■ □ 

From the definition in (j2.12p and Lemma 12.31 the jump conditions of R{z) given in the follow propo- 
sition are easily verified. 

Proposition 2.4. Let Tir he the oriented contour shown in FigureU^ Denote by R+{z) and R^{z), 
respectively, the limiting values of R{z) on S/j taken from the left and from the right of the contour. 
The jump matrix Jr{z) := R-{z)~^ R-i-{z) has the following explicit expressions. For Ke z = I and 
Imz G {—6,6), we have 



Jr{z) 

On the positive real line, we have 



l-A±(z)V±(z) V±(z) 
-A±(z) 1 



Jr^x) 



for X £ (0, 1), and 



Jr{x) 

for x G (l,oo). Furthermore, we have 



1 V- 




(x) - v+(x; 
1 



Jr{z) 



1 

A±(z) 



for z = i Im z with Im z G (—5, 6) and z = Re z ziziS with Re z G (0, 1), and 

Jr{z) 



1 v±(z; 

1 



for z = Re z zizi6 with Re z G (1, oo). 
For simplicity, we define 



e{z) 

viz) 
C 

W{z) 



UTTZ — /3vr/2, 

— zlogc, 

2iTTc-^l''n^, 

(nz) ^-^r(nz + /3/2) 



r(nz + 1-/3/2) 

In view of ()2.4p and the above notations, the functions defined in (j2.13p and (|2.14p become 

Cz^-^WB^ , . 2isin^e±^^+™ 



V± 



2i sin6'e=F»«+^ 



and 



Ah 



(2.15) 



(2.16a) 



(2.16b) 



(2.17a) 



(2.17b) 



(2.18) 
(2.19) 
(2.20) 

(2.21) 



(2.22) 
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id 



Q{z) 



-iS 



Q{z) 


( 1 OA 

^-A+(z) \) 






(I -v+(z)A 
KO 1 ) 




Q{z) 


( 1 OA 

V-A_(z) \) 




1 

Q{z) 


(I -V-{z)\ 
KQ 1 ) 





Q{z) 



Figure 1: The transformation Q ^ R and the oriented contour Tifi.- 
It is easy to see that 
for z £ C±. Also, 
for z = X G (1, oo), and 



A±V± = e^^'' 

4 sin^ 6 



A_ - A, 



for z = X G (0, 1). 



Proposition 2.5. T/ie matrix-valued function R{z) defined in Ii2.12\) is the unique solution of the following 
Riemann-Hilbert problem: 

(Rl) R{z) is analytic in C \ S/j; 

(R2) for z G T,ji, R+{z) = R^{z)Jji{z), where the jump matrix Jr{z) is given in Proposition \2.4\ 
(R3) forzeC\ ^R, R{z) = 1 + 0{\z\-^) as z^oo. 

Proof. It follows readily from Proposition 12.21 and the definition of R{z) in (j2.12p . ■ 

3 Some auxiliary functions 

To construct our parametrix, we should introduce some auxiliary functions. First, define the two constants 

~- - ^ + ^ (3.1) 



a :- 



and b :- 



1 + Vc ' ' 1 - Vc 

These constants are the two turning points for the Meixner polynomials; see [H (2.6)]. Let 



\/bz — 1 + y/az — I ^z — \/ z — b 

[z) := z log - log ^= 

\Jbz — 1 — y az — 1 v-2 — a — v -z — o 



for z G C \ (— oo, 6] and 



, , , a/1 - az + \/\-bz \Jb - z + yja- z 

[z) := z log , , - log ■ 



az 



— bz \Jb — z — ^ a — z 



(3.2a) 



(3.2b) 
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for z e C \ (—00,0] U [a, 00). These two functions are analogues of the ^-function and (/)-function in 
It is clear from the defintions that 



(z) ± i7r(l — z) 



(3.3) 



for z G C±. As 2; —7- 00, we have 



, . , ^/b + ^/a b-a 1, 

(z) = 2;log^= -logz + log— 1 + 0(-, 

vo — ^/a 4 z 



Here we have used the fact that ab = 1. Put 



/:=2W^-2, 



(3.4) 



and recall the definition of v{z) in ()2.19p . Since {Vb+ ^/E)/{Vb - Va) = 1/Vc by ([31]), it follows from 
the above two equations that 



{z) + v{z)/2 + 1/2 = log z + O(-) 



as z — )■ 00. For convenience, we define 



F(z): 

Note by that </>(6) = and 



\n(t){z) 



2/3 



and F{z) :- 



-^n(/)(z) 



2/3 



(3.5) 



(3.6) 



,,, , -, \/bz — 1 + \/az — 1 , — 6z + \/l — 

^ (z) = log , = log ■ 



\/6z — 1 — — 1 — 6z — "v/F 



az 



(3.7) 



Using (3.3), it is readily seen that <j){a) = and 0(0) = ^ logc. The mapping properties of the functions 
(/)(z) and 0(z) are illustrated in Figure [2j From this figure and the definitions ()3.2p and p.6p . we have 
the following proposition. 



Proposition 3.1. For z € C \ [a, 6], we /laue 

±argF(z) G (— 7r,7r) and argF(z) E (— 7r,7r). 
For Re z S (a, 6) and itlmz G [0,(5], we have 

± arg F(z) G (7r/3, vr] and =F arg F(z) G (vr/S, vr]. 

For X > and 5 sufficient small, we have 



(3.8) 



(3.9) 



Re 4>{x lb i(5) ~ < 



—25 arctan 
(l){x) — ttS, 



1 — ax 
bx — 1 



X > b\ 



, a < X <b; 



< x < a. 



(3.10) 
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Proof. It is easy to prove (|3.8p and (|3.9|) by using (j3.2p . (j3.6|) and Figure [21 For small 5 > 0, we have 
from a two-term Taylor expansion 

Re (j){x ± i5) ~ Re <l)± (x) =F (5 Im (x) . 



Note by ((SZD that 



Im (/)'_!_ (x) = < 



0, 

±2 arctan 



1 — ax 
bx-1 



X > b; 
, a < X < 6; 



< X < a. 



Moveover, (|3.2p and (|3.3p imply 



(/)(x), X > 6; 



Re</)-i-(x) = < 



0, 



a < X < b; 



^(p{x), < X < a. 
Thus, (jS.lOp follows from the above three equations. I 
Define 



N{z) :-- 



( (^_l)^( v^+v^ )/3 

{z - a)i/4(2 _ 6)1/4 



n 



(z-a)i/4(2_ 6)1/4 



\ (z - a)i/4(2 _ 6)1/4 (^_a)i/4(^_ 6)1/4 / 

It is easy to verify that N(z) is analytic in C \ [a, b] and 



N+{x) = N_{x) 



-|x-l|''-i 
|x-lp-^ 



X G (a, 6). 



(3.11) 



(3.12) 



The matrix N{z) is analogous to the matrix N{z) in [5l|T5]. Now we introduce the Airy parametrix which 
is also similar to the one in [5| [15] . For z € C± , define 



^(,^ ._ I Ai(z) -iBi(z)\ A Tl/2\ 
^l^i - l^^Ai'(z) Bi'(z) VO 1/2 J 



It is clear that 



A+(x) = A_(x) 



1 -1 
1 



(3.13) 



(3.14) 



on the real line. For convenience, set oj = e^'^*/^. Note that (cf. [1, (10.4.9)]) 

2uj A[{ujz) = - Ai{z) + iBi{z) and 2a;^ Ai(tj^z) = - Ai(z) - i Bi(z) 
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We obtain from (|3.13|) 



A{z) 



' f Ai(z) uj^Ai{uj^z)\ 

^iAi'(z) iuj Ai' {uj^ z) J 

Ai[z) —ujAi{ujz) 

\iAi\z) -iu?A]!{uJz) 



(3.15) 



Furthermore, in view of (cf. [1, (10.4.7)]) 

Ai{z) + uj A\{ijjz) + ij? A\(u?z) = 0, 

we have 



A{z) 



1 

±1 1 



-bjA\(ujz) b? A\{(j?zy\ 

-iu? A{{ujz) iijjA{(b?z)j 

-u? A\{tJ^z) —bjA\(bJz) 

-ibj A{ (ijj^ z) —iijJ^A^(ijJz) 



(3.16) 



Recall the asymptotic expansions of the Airy function and its derivative (cf. [121, p. 392] or [16:^ p. 47]) 



Ai(z) 



-1/4 



2V^ 



-^23/2 

e -i 



Ai'(z) 



=e 3 ^ 



2^ 



s=0 ^3 



3/2-)s 



(3.17) 



as z — )• oo with | argz| < vr, where Us,Vs are constants with uq = vq = 1. For argz G (— 7r,0], we have 
aicg{u}z) G (— 7r/3, 27r/3]. Thus, by using (j3.17p we obtain as z — )• oo with argz G (— 7r,0], 



-a;Ai(a;2;) ~ 
-icj^ Ai'(a;z) 



-iZ 2 3/2 



2^^ 



-63" 



2V^ 

ia;^(a;z)-^/'* ^_2((^2)3/2 z^^^ 2^3/2 



2^^ 



-e 3V 



2V^ 



63- 



For argz G [0,7r), we have arg(6t;^z) G [47r/3, 77r/3). Here, we cannot use (|3.17p with z replaced by uj'^z. 
However, since io'^z = io~^z and arg(a;^^2;) G [— 27r/3, 7r/3), we can use p.l7p with z replaced by ui~^z 
and obtain, as z — )■ 00 with argz G [0, vr). 



CO Ai{uj^ z) = u}'^ Ai{uj z) 



-1/4 



U [OJ Z) I _2/-l^-^S/2 —IZ ' 2^3/2 



2^^ 



-6 3\ 



2\AF 

,1/4 



-63' 



• A-// 2 N • A-// -1 N -iL0{uj ^z)^/^ _2c^-l2-,3/2 Z^'^ 23/2 
Al (W^Z) = Al (W ^z) ^ -=r^ 6 3^'^ ^> ~ —63^ 

2^y^: 2-y/7r 



Applying (|3.17p and the above four formulas to (j3.15p gives 



A{z) 



2V^ 



(/ + 0(|z|-='/^))e-i 



i^«/V3 



(3.18) 
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as z — )• oo with argz G (— vr, vr). For argz G (vr/S, vr], we have arg(zi; z) G (— vr, — vr/S] and arg(tt; z) G 
(— 7r/3, 7r/3]. Thus, by using (j3.17p we obtain as z — )■ oo with argz G (vr/S, vr] 



-a;Ai(a;2;) = — cjAi(a; 



icj Ai (ojz) = —iuj Ai'(a; z) 



2^ 



Z ' _2^3/2 

ia;2(w-2z)V4 _2( 2^)3/2 -iz^^ _2 3/2 



20F 



-e 3V 



2V^ 



-e 3" 



and 



Ai(a;^z) = cj^ Ai(a; ^z) 



ia; Ai'(a; z) = iuj z) 



Uj'^ju ^Z) 1/^ ^,2(^-1^)3/2 



-1/4 



20F 



23/2 

63^= , 



2^^ 



2^^ 



20F 



For argz G [— vr, — 7r/3), we have arg(a;^2;) G [7r/3,7r) and arg(c<j2;) G [— 7r/3, 7r/3). Thus, as z — )• 00 with 
argz G [— vr, — 7r/3), we obtain from ()3.17p that 



and 



-up' a:\{u?z 

-iujK]!{uPz) ^ 

— a;Ai(a;z) ~ 
—iup Ai'(a;z) 



-1/4 



2^ 



iu;(a;^z)i/^ ^_ 2(^2^)3/2 



2^^ 

-zzV4 



_2^3/2 
e 3^ , 



20F 

-a;(tJz)-i/4 



20F 



-^^3/2 

e 3^ , 



2V^ 

iuj'^iujz 



e 3 



2(^^)3/2 ^ -iZ 1/^ ^2^3/2 



2/',.,.)l/4 
2^^ 



2^^ 

1/4 



e 3 



(^,)3/2 ^ ^^,3/2 



2^ 



Applying the last eight formulas to (|3.16p gives 



Mz) 



1 
±1 1 



2-^3/4 

2^ 



-;]{I + 0{\z\~^l'))e-i 



(3.19) 



as z —7- 00 with I arg z\ G ('7r/3, vr]. Here the sign it means plus sign when z is in the upper half plane, and 
minus sign when z is in the lower half plane. Finally, we introduce a crucial function which enables us to 
obtain global asymptotic formulas without any cut in the complex plane; see a statement in the second 
last paragraph of Section 1. For z not on the imaginary line, we define 



''e"^r(nz-/3/2 + l) 



D{z) := { 



27r(nz)"^+(i-/5)/2 
2^(_nz)-"^+(/3-i)/2 



Re z > 0; 



Re z < 0. 



(3.20) 



e-"^r(-nz + /3/2) 
The jump of D{z) on the imaginary line is given by 

Jd(z) := D_{z)-^D+{z) = ^2isin7r(nz - /3/2)e±*^('^^-/^/2) = 1 - e±2*^("^-/3/2), (3.21) 
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where D^{z) (D_{z)) is the hmiting value of D(z) taken from the left (right) of the imaginary axis. It will 
be seen in the proof of Lemma [4.1l that the usage of D{z) is to cancel the jump 1 — e='=2«'f("2-/3/2) g^^j-Qgg ^j^g 
imaginary axis. The explicit formula of D{z) is obtained by solving a one-dimensional Riemann-Hilbert 
problem and calculating a Cauchy integral. As n — )• oo, applying Stirling's formula (cf. [Tl, (6.1.40)]) to 
([320]), we have 

D{z) = l + 0{l/n) (3.22) 
uniformly for z bounded away from the origin. 



4 Construction of parametrix 



For Rez ^ [0, 1] or Im^ ^ [—5, 5], we define 



R(z) := V^[Ce''Y'^^N{z)iz - 1)^"^ Q F{zY^/^ A{F) 



CzP-^B{zY 



(4.1a) 



For Rez G (0, 1) and Imz G (—5, 5), we define 
R{z) , 



(-l)'^+i^/^[C7e'^']'"3/2^(z)(l - z)^"' { ] M ^(z)-'^^/^! A{F)ai 



4sin2 0(z)Z)(z)2e"^(^) 
Czf^-^B{zf 



(4.1b) 



Note that R{z) has jumps across the negative real axis and the imaginary axis; they are caused by the 
functions z^~^ and D{z) respectively. Our parametrix R{z) is analogous to that in (4.54)]. The main 
difference is that the parametrix here has a factor involving the auxiliary function D(z) defined in (j3.20p . 
This factor will make our asymptotic formulas valid in much bigger regions, one of which includes the cut 
(—00,0]. As 2 — )• c«, a combination of (j3.5p . (|3.6p and (|3.18p gives 

V^(^. F(z)'^3/^A(F) = e-"'^"3[j + o(i)] =e(-^-fV3^n-3[j + o(i)]. 



Furthermore, it is easily seen from the definitions (j2.8p . (|3.1ip and ()3.20p that we have respectively 
B{z) ~ z", N{z) ~ / and D{z) ~ 1 as z — )• oo. Thus we obtain from (j4.1ap and the above formula that 



R(z) = [C7e"']'^^/2^V'^3[j + (9(i)]e(-^-fV3^na3 

z 



Cz^-^z^'^ 



= i + o{-] 

Z 



as z — )■ oo. Define 

K{z) := {Ce''^)-''^/^R{z)R{z)-\Ce'^^ f^/'^. 
It is clear from (j4.2p and Proposition 12.51 that 

K(z) = I + 0(-) 



(4.2) 



(4.3) 



(4.4) 
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as z — 7- oo. Let T,^ denote the oriented contour consisting of S/j in Figure 1, the negative real axis, and 
the two infinite hnes from z = zizi6 to z = zizioo on the imaginary axis. The jump matrix of K{z) is given 

by 

Jk{z) := K^{z)-^K+{z) = {Ce^Y''^l'^R^{z)jR{z)R+{z)-^{Ce^^Y^''^ . (4.5) 
Lemma 4.1. Jk{z) = I + 0{l/n) and K{z) = I + 0{l/n) as n ^ oo. 

Proof. In view of the structure of the contour S/^, we divide our discussion into eight cases and consider 
each case separately. 

Case I. For Rez = 1 and ±lmz G [0, 6], we have from (|2.15p and (12.22P that 

/ , o.. Cz^-^B^W \ 



Jb{z) 



1 - e 



±2ie 



- 2i sin 9e^'^+' 
V Czt^-^B^W 
This together with (j4.ip and (j4.5p gives 



2i sin 6'e=F»^+™ 
1 



/3-1 



Jk{z) = N{z){z-1)- 



1 i 



0-3 



1 i 



A F{zy^/^ A{F) 



2 sin (9 



2 sin 6* 



n -1 



Mf(z)-'^3/ViA(F)<ti 



1-/3 



(1 - z)—'''-'N{z)-^{-l) 



(4.6) 



On account of (|3.6p . we obtain from (|3.19p that as n — )• oo, 



1 i 

i 



;jF(.r/*A(F) = [/+o(i)i£^(^\ » 



and 



1 i 



Mi7(^)-3/V^A(F)ai 



1 



From (l2l8D and ([S^D we have 

As n — )• oo, applying Stirling's formula (cf. [1, (6.1.40)]) to ()2.2ip yields 

W{z) = l + 0{l/n) (4.7) 

uniformly for z bounded away from the negative real axis. Applying the last four equations and (|3.22p to 
()4.6p gives 

1 ±(1 - I)2^)g-2n^Ti7r/3 



Jk{z) = N{z){z-1)^ 



=F(1 - D-'^W-^)e^'"l' 1 - e=^2ie(2 - D'^W - D-'^W 



x(z- l)^"3iV(z)-H^ + 0(-)] 

n 



I + Oi- 



n 
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Here we have used the fact that in the present case, Re (j){z) < and Re (p^z) > 0; see Figure [21 
Case II. For z = X £ [1, oo), we have from ([2J6]l and (fr22]) 



1 -Cx^-i^^p^e"™ 
1 



This together with ([3H]) . (fiJ]) and ([33]) gives 



Ji4-(x) 



iV(x)(x - 1) 



. / 1 i 



1 1-Z)2VF 
1 



N{x)ix - 1] 



/3-i^„ / I i 



-(73 



(4.8) 



Note that the matrices N(z) and F{z)"'^^^ are both discontinuous across the interval [1, b). But a combi- 
nation of them makes the jumps vanish. Observe from Figure [2] and (13. 6p that aTg(j)±{x) = ib37r/2 and 
argF±(a;) = ibvr for x G [1,6). Thus, we have F+(x)'^3/4 ^ F_ (x)'^3/4gi7ra3/2^ ^^^^^ follows from (fXT2]) 

that the matrix 



Niz)iz-1) 



-0"3 



[)F{zr" 



has no jump on the interval [1,6). Applying (|3.22p and (jiT]) to gives Ji<:(2;) = / + 0(l/n). 

Case III. For z = x G [0, 1], we can proceed in a similar manner as in Case II and obtain Jk{x) 
l + 0{l/n). 

Case IV. For z = Rez±i5 with Rez € (l,oo), we have from (2.17b) and (|2.22p 



Jr{z) 



This together with (14. ip and (14. 5p gives 



Ji^(z) = iV(z)(z-l)- 



1 z 



2i sin (9e=F«^'+" 



1 i 



* J F(z)'^^/4 A(F)e"<^'"3 



2i sin 6*6^*^ 

Vo 1 



M F(z)'^3/^A(F)e"'^'^=' 



-1 



(z- l)V'"3iV(z)-^ 



(4.9) 



Note from p.Sp that in this case, we have argF(2;) € (— 7r,7r). Thus, coupling (j3.6p and (j3.18p . we obtain 

1 „,1. 



1 i 

i 



: F(z)'^^/^A(F)e"^"3 = r/ + o(_)]. 
1/ n 



Applying this to ()4.9p yields 



Jk{z)=I + 0{-). 

n 
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Here we have used the facts that =p2i sin 6e^^^ = 1 — e^^'^ ~ 1 as n — )• oo and Re(— it i9) = — n(Re (j) + 
7r6) < 0; see ([2J8]l and (fXTO]) . 

Case V. For z = x it i(5 with x G (0, 1), we have from (2.17a) and ()2.22p 



jRiz) 



This together with (|4.ip and ()4.5p gives 



Jx(^) = (-ir+iiv(z)(i-z)^-« 



2i sin 6*6=^*^+"'' 



1 i 

i 1 



2 sin 6* 

iD-2^-ie=t*^ (2sin6')-i 



1 i 
i 1 



1-/9 



(4.10) 



Note from (|3.8p that in this case, argF(z) G (— 7r,7r) and argF(z) G (— 7r,-7r). Thus, we obtain from ()3.6 
and (IXTHD 



1 i 

i 



! ) F(z)^'^3/Vi A(F)ai = [/ + 0(-)]^^ 
1/ n 



and 



1 i 



H -1 



n 



^ ljF{zr/'A{F) 

Applying the last two equations to (|4.10p yields 

/ 2(— 1)"+^ sin^e"'^""'^'^*'^*^"^^^''/^ 



1 



/9-1 



Jx(z) = N{z){z-l)- 



0-3 



;(_1)"+1 2)^2p^-lgn0+n0±«e=F«7r(l-/3)/2 



\ 

gn(/>— n(/>=Fi7r ( 1 — /3) /2 

2(-l)"+isin6l / 



x(z- l)V^'^3^(z)-i[/ + 0(i)]. 

n 

Using (2.18), (3.2), (3.3) and (3.10), one can show that 

Re{n0 + n^±ie} = Re{2n(^} < 

and 



2(-l)"+^ sij^0gn<A-n</,±m(i-/3)/2 ^ sin Se^^^ = 1 - e±2if? ~ 1. 



Thus, we again have 



Jfc(z) = / + 0(-), 
n 



as n — )• 00. 
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Case VI. For z = ±iy with y G (0, 6), we have from pTT]) and (fX^ 

/ 1 0\ 



Jr{z) 



This together with (|4.ip and (|4.5p gives 



Jx(^) = (-ir+^iv(z)(i-z)^'^^ 



1 i 



2i sin 6»e=t*^+"'' 



/2i:»_ sin( 



1 i 



M i7(2)'^3/4A(F) 



MF(z)-'^3/ViA(F)ai 







\ 2D_sin6l/ 



(4.11) 



Note from ()3.8p that in this case, argF(2) G (— 7r,7r) and argF(2) G (— 7r,7r). Thus, as in Case IV we 
have from ([MD and (fXTH]) 



1 i 

i 



■I \ ~ u ~ 1 e'^'l"^^ 



and 



1 i 



^ ;jF(z)'^3/4A(F) 

Applying the above two equations to (|4.1ip yields 

/ 2(-l)'^+isin 



n 



Jk{z) = N{z){z-\)^''-'' 



Jj^^n(t>-n4>TiTT{l-l3)/2 



0{e 



2n Ret 



Jj^Qn<p-n<t,TiTr{l-l3)/2 

2(-l)"+isin0 y 



x(z- l)^'^«Af(z)-i[/ + 0(i)]. 



Here we have used (2.21), (3.20) and the asymptotic formula for r(x it iy) as y — )• +cxd. Since 

2(-l)'^+i sin^e'"^-''^^*^^^-^)/^ 



and Re(/)(2;) < 0, as before we again have Jk{z) = I + 0{l/n) as n — )• oo. As mentioned in a statement 
following (3.21), the usage of D{z) defined in (15:20]) is to cancel the jump 1 - e^'^'^ = 1 - e±2i'r(nj;-/3/2) _ 
Without this function, the jump matrix Jk{z) is not asymptotically equal to the identity matrix in this 
case. 
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Case VII. For Rez = and | lmz\ > 6, we have Jr{z) = /; see Figure 1. Thus, (j4.ip and (|4.5p imply 

Jk{z) = 



N{z){z - 1) V'^a f 1 ^ j F(z)'^3/4 



T-o-3 



(4.12) 



Note that by (j3.2ip . J/) = 1 — e±2i7r(n2-/3/2) exponentially small for | Imz| > 6. From (|4.12p . it again 
follows that Jk = I + 0{l/n) in this case. 

Case VIII. For z = x G (— oo,0), we have Jr{x) = I. Thus, (|4.ip and (|4.5p imply 



0-1 i7r{l-/3) 

Jj^{z) = N{x){l- x)—''^e^^''^ 



1 i 



1 i 



„i7r(/3- 1)0-3 



1-/3 i7r(l-/3) , 

(1 - x)— '"»e^^'"3iV(x)-^ 



(4.13) 



Here we have used the fact that e''(<^+-'^-) = e^™'^ = 1 for x < 0; see Note from ([3:8]) that 

argF-t S (— vr, vr) in this case. Hence, we obtain from (j3.6p and (j3.18p 

] F±(x)'^3/4 A(i7^)^n0±.3 = 1 [J + 0(i)]. 



Applying the last equation to (I4.13|) yields 

13-1 i7r(l-i9) . ,x !7r{l-/3) 1- fi -, 1 1 

J^(z) = Ar(x)(l-x)V-3e^T^-3gM/3-i)-3g^T^-3(i_^)^-3^(^)-i[j + C>(_)] = j + 

In conclusion, we have shown that Jk{z) = I + 0{l/n) on the contour of K{z). It is not difficult 
to verify that the multiplicative cyclic condition (3.30) in [14] holds for the jump matrix Jk{z). An 
application of Theorem 3.8 in [14] then gives K{z) = I + 0{l/n) as n — t- oo. ■ 

Combining Lemma O with ([23]), (USD, ^M), (|2l^ and (g^l), we obtain 



7r„(nz - 13/2) = vJ'B{z)Rii{z)[I + 0(-)] 



n 



for Rez ^ [0, 1] or Imz ^ [—5, 5], and 



7r„(nz - /3/2) = n^B{z)[Rii{z) + A±(z)iii2(z)] [/ + 0(-)] 

n 



(4.14a) 



(4.14b) 



for Rez G (0,1) and Imz G (0,±(5). 



5 Main results 

Theorem 5.1. As n ^ oo, we have 



7r„(nz-/3/2) 



n"'^/^^D{z}e 
( 



rfu(2)/2+n//2 



l+V 2 — 6 



r + 



z—a—Jz—h \ 



zW-^)/2(z - a)i/4(2 _ 6)i/4F(z)-V4 



Ai(F) 



^(/3-l)/2(^ _ a)l/4(^ _ 5)1/4^7(^)1/4 Al'(i^)| [1 + 0(^)] 



(5.1) 
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for Re z ^ [0, 1] or Im z ^ [—5, 5], and 

7r„(nz-/3/2) = (-n)"V7?L>(z)e"^(^)/2+n//2 

^ \/b~z+^ a—z _|_ ^ \/b—z—^/a—z 



^{/3-l)/2(5 _ ^)l/4(o _ 2)1/4 i7(^)-l/4 



[cos6iAi(F) -sm6'Bi(F)] 



(5.2) 



+ 



if 



zW-i)/2^b - z)V4(a _ z)l/4i?(2)l/4 



[cos 9 Ai'(F) - sin 6 Bi'(F)] I [/ + 0(-)] 



n 



/or Rez E (0,1) a^^^i Imz S {—6,6), where the constant I is given in (3.4) and the functions v{z),F(z) 
and D{z) are respectively given in (2.19), (3.6) and (3.20). The asymptotic formula on the boundary of 
the two regions can he obtained by taking limits from either side. 



Proof. From (IXTTD . (IXT3D . (liH) and (I4.14ap . it is easy to obtain ([5l]). We now prove ([52]) • Define 



Q{z) := R{z) 



1 0\ 
A±(z) 1 ■ 



Prom (I4.14bp . we have 



TTninz - (3/2) = n^B{z)Qu{z)[I + 0{-)]. 

n 



(5.3) 



(5.4) 



Thus, we only need to calculate Qii{z). First, we observe from (j2.22p . (j4.ip and (|5.3p that 



Qiz) 



{-l)''+^V^[Ce"Y^/^N{z){l - z] 2 
/ 1 0\ 

j^^±id{z) 

1 



- / 1 i 



i 1 



F(z)~'^^/Vi A(F)cJi 



\2D'^{z)W{z)s\ne{z) J 



4sin2 6'(z)L>(z)2e™(^) 
C^/3-iB(z)2 



fT3/2 



Second, (jS.lip gives 



iV(z)(l - 2)^-^3 



1 i 
i 1 



/ ( V5^+v^ )/3 _ (. 



/b—z—^/a^ 



If 



Ib—z+y/a^ 



[)/3 + ( v^-v^ )/3 \ 



z)l/4(a _ 2)1/4 



-i(6-z)i/4(a-z)i/4 



-i(6-z)i/4(a_ 2)1/4 



(6-z)i/4(a_2)i/4 
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Finally, (|3.13|) implies 



CTl A(F)f7i 



\2D'^W 



1 



0\ 



1 



(2 sin 6*) 



sm ( 



Bi'(F2 iA^{F)\ (1/2 
-iBi(F) Ai(F) 



2 sin (9 







±1/2 ly \iD-'^W-^e^''^ (2sin6l)-i 

iAi'(F)\ 
2 sin 61 



sineBi'(F) - [cos^ + _ l)e±*^] Ai'(F) 



o 1 _L-^ - Ai(F) 

\ -i{sin 6* Bi(F) - [cos 6 + (D-^^^-i - l)e^'^] Ai(F)} — -- / 

\ 2sin0 / 

Applying the last three equations to (|5.4p gives 



7r„(nz-/3/2) 



nt)(2)/2+ri«/2 



/6— 2+v'<^^2: \/3 

2 



2 ) 



^(/3-l)/2(5 _ ^)l/4(^ _ ^)l/4^(^)-l/4 

X [cos 61 Ai(F) - sin 9 Bi(F) + e^''%D~^W~^ - 1) Ai(F)] 



2 



^(/3-l)/2(5 _ ^)1/4(q _ 2)1/4^7(^)1/4 

X [cos e Ai'(F) - sin 9 Bi'(F) + e^'^fD^^^-^ - 1) Ai'(F)l \[I + 0(-)]. 

J n 



On account of (j3.2p . (j3.6p . and (jS.lTp . Ai(F) and Ai'(F) are exponentially small when z approaches the 
origin; by (|3.22p and (j4.7p . we also have D~^W~^ — 1 = 0(l/n) for z ^ 0. Since we can always neglect 
the terms {D-'^W-^ - 1) Ai(F) and {D-'^W-'^ - 1) Ai'(F), formula (5.2) is proved. 

To justify that the asymptotic formula on the curve separating the two regions can be obtain by 
taking limits from either side, we just note that the regions of validity of both formulas (5.1) and (5.2) 
can be slightly extended beyond their boundaries, and that in the overlapping region these two formulas 
are asymptotically equal. ■ 

Remark 5.2. We would like to mention that our results coincide with those obtained in The 
formulas (6.9) in [8] and (2.35) in [9] are asymptotically equal to l[5.1\) in the present paper, while the 
formulas (6.27) in [8] and (4-19) in [9] are asymptotically equal to ^5.2\) . 
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